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Abstract 

Observations using the Rossi X-ray Timing Explorer (RXTE) have discovered 
dozens of accreting neutron stars with millisecond spin periods in low-mass bi- 
nary star systems. Eighteen are millisecond X-ray pulsars powered by accretion or 
nuclear burning or both. These stars have magnetic fields strong enough for them 
to become millisecond rotation-powered (radio) pulsars when accretion ceases. Few, 
if any, accretion- or rotation-powered pulsars have spin rates higher than 750 Hz. 
There is strong evidence that the spin-up of some accreting neutron stars is limited 
by magnetic spin equilibrium whereas the spin-up of others is halted when accre- 
tion ends. Further study will show whether the spin rates of some accretion- or 
rotation-powered pulsars are or were limited by emission of gravitational radiation. 
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1 Observed Properties of Millisecond X-ray and Radio Pulsars 



In this brief review, a pulsar is considered a millisecond pulsar (MSP) if its 
spin period P s is < 10 ms (spin frequency v s > 100 Hz). Seven accretion- 
powered and 13 nuclear-powered X-ray MSPs are now known (see Table 1). 
All are accreting neutron stars in low- mass X-ray binary systems (LMXBs). 
Two of the accretion-powered MSPs become nuclear-powered MSPs during 
thermonuclear X-ray bursts. Kilohertz quasi-periodic oscillations (QPOs) have 

* This research was supported in part by NASA grant NAG 5-12030, NSF grant 
AST 0098399, and the funds of the Fortner Endowed Chair at the University of 
Illinois. 



Preprint submitted to Elsevier Science 



5 February 2008 



Table 1 
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"Spin frequency inferred from periodic or nearly periodic X-ray oscillations. 
A: accretion-powered millisecond pulsar. N: nuclear-powered millisecond pulsar. K: 
kilohertz QPO source. See text for details. 



so far been detected in the accretion-powered X-ray emission of more than 
two dozen neutron stars in LMXBs, with frequencies ranging from ~ 100 Hz 
to ~1300 Hz (Lamb 2003; van der Klis 2006). They have been seen in 8 of 
the 13 known nuclear-powered MSPs and 2 of the 7 known accretion-powered 
MSPs (see Table 1). 

The discovery of a pair of kilohertz QPOs in the 401 Hz accretion- and nuclear- 
powered X-ray pulsar SAX J1808.4-3658 (Chakrabarty et al. 2003; Wijnands 
et al. 2003), of a pair of kilohertz QPOs in the 191 Hz accretion-powered X- 
ray pulsar XTE J1807.4— 294 (Linares et al. 2005), and of burst oscillations 
in the 314 Hz accretion-powered X-ray pulsar XTE J1814— 338 (Strohmayer 
et al. 2003) has confirmed the relationship between these three distinct types 
of high-frequency X-ray oscillations. 
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Except during the first few seconds of some bursts, the nuclear-powered oscil- 
lations of SAX J1808 and XTE J1814 have the same frequency and phase as 
the accretion-powered oscillations. These observations establish beyond any 
doubt that (1) the nuclear- and accretion-powered oscillations are both pro- 
duced by spin modulation of the X-ray flux from the stellar surface and (2) the 
magnetic fields of these stars are strong enough to channel the accretion flow 
and enforce corotation of gas heated by nuclear burning, which requires fields 
> 10 7 G (Miller, Lamb & Psaltis 1998; Lamb & Yu 2005). Conversely, the 
nearly sinusoidal waveforms and low amplitudes of these oscillations indicate 
that the surface magnetic fields of these neutron stars are less than ~ 10 10 G 
(Psaltis & Chakrabarty 1999). These results confirm that the frequencies of 
both the accretion- and the nuclear-powered X-ray oscillations are the spin 
frequencies of these neutron stars. 

None of the 18 currently known accretion- and nuclear-powered X-ray pulsars 
have spin frequencies higher than 619 Hz (see Table 1). There is no known 
observational bias against detecting such pulsars. This indicates that few, if 
any, such pulsars have spin frequencies higher than 750 Hz (Chakrabarty et 
al. 2003; M. C. Miller 2005, personal communication). Although too few such 
pulsars have been observed to establish their spin distribution accurately, the 
spins of the known MSPs nevertheless provide valuable information about 
their production and evolution, as discussed below. 

A weak-field accreting neutron star typically produces a pair of strong kilo- 
hertz QPOs, with frequencies that vary by as much as a factor of five as the 
star's X-ray flux varies by a factor of two or three (see Miller et al. 1998; Lamb 
2003; Lamb 2005; Lamb & Miller 2006; van der Klis 2006). As the frequencies 
of the two kilohertz QPOs vary by hundreds of Hertz, their separation Az/qpo 
remains roughly constant. All the kilohertz QPO pairs detected in X-ray pul- 
sars have Az/qpo approximately (sometimes very accurately) equal to either 
v s or 0.5 v s . For example, the kilohertz QPO pair discovered in XTE J1807.4 
has Az/qpo — v s (Linares et al. 2005) whereas that discovered in SAX J1808.4 
has Az/qpo — 0.5 v s (Wijnands et al. 2003). These findings show that the spin 
of the star plays a central role in generating its kilohertz QPO pair. This is 
possible only if the neutron stars that produce kilohertz QPOs have magnetic 
fields > 10 8 G (Miller et al. 1998; Lamb & Miller 2001; Lamb & Miller 2006). 
The central role played by the star's spin frequency means that information 
about its spin frequency can be extracted from the frequencies of its kilohertz 
QPOs, even if accretion- or nuclear-powered oscillations have not yet been 
detected. 

Altogether, more than two dozen accreting neutron stars have spin rates 
and magnetic fields high enough for them to become rotation-powered radio- 
emitting MSPs when accretion ceases, supporting the hypothesis (Alpar et 
al. 1982; Radhakrishnan & Srinivasan 1982) that such neutron stars are the 
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progenitors of the radio MSPs. Many new rotation-powered radio MSPs have 
been discovered in recent years and dozens are now known (see Lorrimer 2005). 
Eighty binary and millisecond radio pulsars have been associated with the 
Galactic disk. More than one hundred radio pulsars have been discovered in 
24 globular clusters in the Galaxy (Lorrimer 2005). In November 2005, a radio 
pulsar was discovered in Terzan 5 with a spin frequency of 716 Hz (Hessels 
et al. 2006), higher than the 642 Hz spin frequency of the first radio MSP 
discovered (Backer at al. 1982). 



2 Production of Millisecond X-ray and Radio Pulsars 

Production and spin evolution of accretion-powered millisecond X-ray pulsars.— 
The timescale on which accretion from a disk doubles the spin rate of a slowly 
rotating magnetic neutron star is (Ghosh & Lamb 1979b; Ghosh & Lamb 
1992) 

/ /iV \ ~ 1+a / 3 

*. , 2^/ [ M (G M rm) v V 10V yy yy . « 

Here v s , M, and / are the star's spin frequency, mass, and moment of inertia, 
M is the mass accretion rate onto the star (not the mass transfer rate), Me 
is the mass accretion rate that produces the Eddington luminosity, and r m is 
the angular momentum coupling radius. In the final expression on the right- 
hand side of equation (1), the weak dependence of t s on M, /, and the star's 
magnetic field has been neglected; a is 0.23 if the inner disk is radiation- 
pressure-dominated (RPD) or 0.38 if it is gas-pressure-dominated (GPD). 

The first papers to argue that pulsing X-ray stars are accreting magnetic neu- 
tron stars also argued that accretion of angular momentum would spin such 
stars up until their angular frequencies are approximately equal to the or- 
bital angular frequency at the disk-magnetosphere coupling radius (see Lamb, 
Pethick & Pines 1973). Accretion from a disk drives the neutron star's spin rate 
upward or downward on the timescale (1) until it reaches the "magnetic spin 
equilibrium" period P eq . Subsequent two-dimensional semi- analytical MHD 
calculations (Ghosh & Lamb 1979a; Ghosh & Lamb 1979b; Ghosh & Lamb 
1992) supported the conclusion that the torque on a magnetic neutron star 
accreting from a disk vanishes at a certain period and provided a more accu- 
rate value of this period. Recent two- and three-dimensional numerical MHD 
calculations (Romanova et al. 2002) have provided further support. 

By 1982, spin-up by accretion in binary systems was being considered to 
explain the short-period (~ 60-100 ms) rotation-powered radio pulsars that 
had been discovered (see Taylor & Stinebring 1986). That year Alpar et al. 
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(1982) and Radhakrishnan & Srinivasan (1982) argued that some neutron 
stars have high enough accretion rates and weak enough magnetic fields that 
they can be spun up to periods as short as a millisecond. That same year 
Backer et al. (1982) discovered a rotation-powered pulsar with a 1.56 ms pe- 
riod (PSR 1937+21). 

Magnetic neutron stars formed in close binary systems are thought to be spun 
down initially, by emission of particles and electromagnetic radiation, and then 
spun up by accreting angular momentum when mass transfer occurs. This 
"recycling" process produces weak-field neutron stars with high spin rates. In 
order to obtain a definite equilibrium spin period ("spin-up line"), Alpar et al. 
and Radhakrishnan & Srinivasan assumed that the accretion rate at the end 
of mass transfer is the Eddington rate, even for low-mass binary systems. This 
is the highest accretion rate that is possible, but it is not the accretion rate 
expected at the end of mass transfer in the relevant binary systems. Indeed, 
there is strong evidence that the accretion rate that produces most recycled 
pulsars is less than the Eddington rate (see below). 

As noted above, few if any accreting neutron stars in LMXBs have spin rates 
higher than 750 Hz. Possible explanations for this include (see Lamb & Yu 
2005) (a) the finite amount of angular momentum that can be accreted during 
binary mass transfer; (b) the limiting effect of the magnetic braking compo- 
nent of the accretion torque, which becomes dominant at high spin rates; and 
(c) the possible existence of one or more other braking effects — such as grav- 
itational radiation — at high spin rates. I discuss these possibilities in turn in 
the following paragraphs. 

(a) Some neutron stars in LMXBs may continue to spin up throughout the 
mass transfer phase of binary evolution but never reach a spin rate higher than 
a few hundred Hz simply because the spin- up timescale given by equation (1) 
becomes so long that mass transfer ends before such a high spin rate can be 
reached (Lamb & Yu 2005). This appears to have been overlooked in many 
analyses. 

A few neutron stars in LMXBs are accreting steadily at rates ~ M E , but most 
appear to be accreting at average rates < 10~ 2 M E . A star spinning at 600 Hz 
with an average accretion rate ~ 10 -3 M£i has a spin-up timescale ~ 2 Gyr, 
even if the magnetic braking component of the accretion torque and emission 
of gravitational radiation are both negligible. This is longer than the mass 
transfer phase of many binary systems. In such systems, mass transfer will 
end before the star can be spun up further. 

(b) The spin rates of many neutron stars in LMXBs appear to be limited 
by the magnetic braking component of the accretion torque (Ghosh & Lamb 
1979a; Ghosh & Lamb 1979b; Ghosh & Lamb 1992; Lamb & Yu 2005). The 
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Fig. 1. Pulsar evolution diagram (PED) used to analyze the evolution of accre- 
tion-powered X-ray pulsars and the production and evolution of rotation-powered 
radio pulsars. The data points are known rotation-powered pulsars; those of pulsars 
in binary systems are encircled. The data are from Hobbs & Manchester (2004). 

properties of accreting neutron stars in LMXBs provide some evidence for this, 
but the most compelling evidence for magnetic spin equilibrium is provided by 
the properties of the rotation-powered radio MSPs that are discussed below. 



The spin and magnetic field evolution of accretion-powered X-ray pulsars, the 
production of rotation-powered radio pulsars, and the spin and magnetic field 
evolution of the latter can be analyzed using the pulsar-evolution diagram 
(PED) shown in figure 1 (Ghosh & Lamb 1992; Lamb & Yu 2005). 

The lines in the PED that slope downward to the right show the magnetic 
dipole braking relation = 3.2 x 10 19 (P S P S ) 1 / 2 G for a rotation-powered 
pulsar with a constant magnetic field of the strength given by the labels of 
the lines. The heavy dashed line that slopes upward to the right shows where 
stars with a magnetic dipole spin-down timescale equal to 15 Gy would lie. 

The lines in the PED that slope upward to the right show the equilibrium 
spin period P eq specified on the horizontal axis for an accretion-powered pulsar 
that is accreting at the rate given by the labels of these lines and has a dipole 
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magnetic field that corresponds to the magnetic dipole braking rate specified 
on the vertical axis (the actual P s of an accretion-powered pulsar that is in 
magnetic spin equilibrium is of course zero). 

The P eq lines in the PED for M = M E and M = 0.1M E have jumps where 
the structure of the disk changes from RPD (lower left) to GPD (upper right). 
In reality the transition is smooth. The lines for M < 0.01M E do not have 
jumps, because for these accretion rates the disk is GPD at r m , even if the 
star's magnetic field is < 3 x 10 7 G. The effects of the stellar surface and the 
innermost stable circular orbit are not shown. 

The spin rate of an accreting neutron star reflects the accretion torque acting 
on it, averaged over the spin relaxation time t s given by equation (1). As noted 
above, a few neutron stars in LMXBs are accreting steadily at rates ~ M E , but 
most are accreting at rates < 10~ 2 M E and consequently have spin relaxation 
timescales > 10 8 years. Determining the average accretion torque over such a 
long interval is difficult, because of our lack of knowledge of the accretion rate 
and magnetic field histories of these stars (Miller et al. 1998; Lamb & Yu 2005). 
Many accrete only episodically. For example, the recently-discovered accretion- 
powered MSPs in LMXBs usually have accretion rates < 10~ A M E , but also 
have outbursts every few years. During these outbursts, which last for several 
weeks, their accretion rates rise to ~ 10~ 2 M E before falling back to < 10~ 4 M E 
(Chakrabarty et al. 2003; Strohmayer et al. 2003). Also, the magnetic fields of 
accreting neutron stars appear to decrease by factors ~ 10 2 -10 3 as a result of 
accretion, perhaps on timescales as short as hundreds of years (see Lamb & Yu 
2005). All these variations make it difficult to determine the average accretion 
torque over long time intervals. If accretion torque theory can be validated, 
the argument can be reversed and spin-evolution measurements can be used 
to determine the average accretion rate over binary evolution timescales. 

The spin frequencies and accretion rates of the 18 known accretion- and 
nuclear-powered MSPs in LMXBs are consistent with long-term magnetic spin 
equilibrium if their dipole fields are between 3 x 10 7 G and 3 x 10 8 G and their 
time-averaged accretion rates are between 3 x 10~ 3 M E and M E . The spin 
rates and visible X-ray oscillations of the accretion-powered MSPs are under- 
standable if they have magnetic fields ~3 x 10 8 G and have been spun up to 
magnetic spin equilibrium by accretion at rates ~ 10~ 2 M E . The higher spin 
rates of many of the nuclear-powered MSPs are understandable if their mag- 
netic fields are less than ~ 10 s G and their time-averaged accretion rates are 
~ 10~ 2 M E . 

If many accreting MSPs are near magnetic spin equilibrium, the ones with 
lower spin rates should have higher magnetic fields and hence stronger X-ray 
oscillations, other things being equal, making them easier to detect. Indeed, 
six of the MSPs in which accretion-powered oscillations have been detected 
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do have spin rates < 500 Hz whereas only one has a spin rate > 500 Hz (see 
Table 1), consistent with this expectation. More MSPs are needed to draw 
firm conclusions. 

The absence of accretion-powered millisecond X-ray pulsars with frequen- 
cies > 620 Hz is consistent with their being in time-averaged magnetic spin 
equilibrium, if their magnetic fields are greater than ~ 3 x 10 7 G and their 
time-averaged accretion rates are less than ~ 10~ 3 M£. These fields and rates 
are consistent with the other observed properties of neutron stars in LMXBs 
(Miller et al. 1998; Psaltis & Chakrabarty 1999; Chakrabarty et al. 2003). 

(c) The spin-up of neutron stars in LMXBs by accretion may be limited by 
gravitational wave emission, if they have weak enough magnetic fields and 
large enough time- dependent mass- or mass-current quadrupoles (Ushomirsky, 
Cutler & Bildsten 2000; Lindblom & Owen 2002). However, at present there 
is no unambiguous evidence for this. 

We close this discussion of accretion-powered millisecond X-ray pulsars by 
emphasizing again that the paucity of accreting neutron stars with spin rates 
higher than 620 Hz may be due simply to the very long time required for 
accretion to produce such high spin rates (Lamb & Yu 2005): for v s = 750 Hz 
and a time-averaged accretion rate ~ 10~ 3 M£, the spin-up timescale is ~ 2.5 
Gyr, longer than the mass transfer phase in low-mass systems. 

Production and spin evolution of rotation-powered millisecond radio pulsars.— 
The initial spin periods of rotation-powered millisecond pulsars produced by 
accretion in binary systems reflect the appropriately time-averaged magnetic 
fields and accretion rates of the progenitor neutron stars toward the end of the 
mass transfer phase (Ghosh & Lamb 1979b; Ghosh & Lamb 1992; Lamb & Yu 
2005). Comparison of the P s -P s distribution of known rotation-powered radio 
pulsars with the pulsar evolution diagram (PED) shown in Figure 1 therefore 
provides important information about how they are produced: 

(1) The current positions of recycled pulsars in the PED are necessarily to 
the right of their positions when accretion ended, because after accretion ends 
they will be spun down by the magnetic braking torque, even if there is no 
other braking torque. If the braking torque on these pulsars is small enough, 
their positions in the PED will not have changed significantly since accretion 
ended and their current positions can be taken as indicating their positions 
at the end of the mass transfer phase. This is the case for pulsars near and 
below the 15 Gy heavy dashed line. 

(2) If the accretion-powered X-ray MSPs that are the progenitors of the 
rotation-powered radio MSPs are in magnetic spin-equilibrium at the time 
accretion ends, the standard theory of disk accretion by magnetic neutron 
stars predicts that there should be no rotation-powered MSPs to the left of 
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the magnetic spin-equilibrium line for M = Me, because (i) this line defines 
the shortest spin period that can be reached by accretion and (ii) any mag- 
netic braking that occurs after accretion ends will drive the pulsar to the right 
in the PED. We emphasize that an RPD model of the inner disk must be used 
for neutron stars with such weak magnetic fields accreting at such high rates. 

(3) The existence of an upper edge to the distribution of rotation-powered 
MSPs in the PED indicates that the accretion torque vanishes at high spin 
rates, i.e., there is an equilibrium spin rate. The agreement of this edge with 
the magnetic spin-equilibrium line predicted for M = M E by the standard 
model of disk accretion by weakly magnetic MSPs (Ghosh & Lamb 1979b; 
Ghosh & Lamb 1992) indicates that (i) the standard model is roughly correct 
for stars with magnetic fields ~ 3 x 10 8 G to ~ 3 x 10 10 G accreting at 
M fa Me, (ii) some of the progenitors of the rotation-powered MSPs had 
accretion rates ~ M E for long enough to reach the magnetic spin equilibrium 
line for M = Me by the time mass transfer ended, (iii) accretion by these stars 
ended abruptly enough that they did not spin down appreciably after reaching 
this spin equilibrium line, and (iv) the spin rates of an appreciable fraction of 
accretion-powered MSPs are limited primarily by magnetic braking. 

The reported P's and P's of two rotation-powered MSPs recently discovered in 
globular clusters (B1821-24 and B1820-30A; Hobbs et al. 2004) place them 
above the magnetic spin-equilibrium line predicted for M = Me- Either the 
intrinsic P s 's of these pulsars are less than shown or the standard RPD model 
of the inner disk does not accurately describe the accretion flows that spun 
up these stars. 

(4) The MSPs that appear near the upper edge of the MSP distribution in 
the PED must have accreted at near-Eddington rates for at least ~ 1 Myr and 
therefore must have accreted at least ~0.1M Q . 

(5) Many pulsars that are well below the magnetic spin-equilibrium line for 
accretion rates ~ 10~ 2 Me have very long magnetic braking spin-down times. 
This implies either (i) they never reached spin equilibrium or (ii) they accreted 
at rates <C M E for a long time before the end of the mass transfer phase. 

(6) The magnetic spin-equilibrium hypothesis predicts that MSPs should be 
rare or absent below the M = 10~ 4 M E spin-equilibrium line in the PED, 
because stars accreting at such low rates generally will not achieve spin equi- 
librium during their accretion phase. The observed P s -P s distribution is con- 
sistent with this prediction. 

(7) The MSPs near the 15 Gyr spin-down line in the PED were produced in 
situ by accretion at rates < 3 x 10~ 3 M£ rather than by spin-up to shorter 
periods by accretion at rates greater than 3 x 10~ 3 M E followed by magnetic 
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dipole braking, because braking would take too long. This result accords with 
the expectation noted above that most neutron stars in LMXBs accrete at 
rates <C M E toward the end of their accretion phase. 

In summary, the properties of the accretion- and nuclear-powered X-ray MSPs 
discovered in LMXBs using RXTE support the idea that such stars are the 
progenitors of the observed rotation-powered radio-emitting MSPs. After be- 
ing spun down to longer periods by rotational energy loss, the X-ray MSPs 
discovered using RXTE were spun up to millisecond periods by accretion, 
eventually becoming the nuclear- and accretion-powered MSPs we see today. 
Some of these may stop accreting before reaching magnetic spin equilibrium or 
may be spun down as accretion ends, but we expect many of them to become 
rotation-powered radio MSPs. 

I thank L. Bildsten, D. Chakrabarty, P. Kaaret, M. van der Klis, M. C. Miller, 
D. Psaltis, J. Swank, and W. Yu for helpful discussions and W. Yu for help in 
preparing the references. 
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